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Abstract 

O ■ In this paper, we study an adaptive finite element method for multiple eigenvalue problems. 

I We obtain both convergence rate and quasi-optimal complexity of the adaptive finite element 

eigenvalue approximation, without any additional assumption to those required in the adaptive 
finite element analysis for the boundary value problem. Our analysis is based on a certain rela- 
tionship between the finite element eigenvalue approximation and the associated finite element 
• boundary value approximation and a crucial property of eigenspace approximation which are 

also presented in the paper. 
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^ 1 Introduction 

> 

\Q I Adaptive finite element computation is one of efficient approach.es to solve partial differential equa- 

. tions and tias been successfully used in scientific and engineering computing. Its numerical analysis 

OO I has been also derived much attention from the mathematical community. Since Babuska and Vo- 

gelius [3] gave an analysis of an adaptive finite element method (AFEM) for linear symmetric elliptic 
I problems in one dimension, there has been much investigation on the convergence and complexity 

of AFEMs in literature (see, e.g., [Q H [H [El [ISl EH |32] and the references cited therein). In the 
context of the finite element approximations of eigenvalue problems, in particular, we note that 
there are a number of works concerning a posteriori error estimates [5j [T4j [211 ESI EH [33] , AFEM 
convergence [TH [13 [HI [HI [23 and complexity [HI [TSl [SD] . Except for the convergence analysis in 
[16], to our best knowledge, there is no any work about convergence rate and complexity of AFEM 
' for multiple eigenvalue problems. The purpose of this paper is to fill in the gap. 

5^ . We understand that multiple eigenvalue problems are topic in science and engineering, such as 

Hartree-Fock equation and Kohn-Sham equation used to model ground state electronic structures of 
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molecular systems in quantum chemistry and materials science, in which hundreds of thousands of 
eigenvalues and their corresponding eigenfunctions are needed, and among these eigenvalues, most 
are multiple [121 US [IHl 1301 HI] • While the central computation in solving Hartree-Fock equation or 
Kohn-Sham equation is the repeated solution of linear Schrodinger type equation, of which adaptive 
finite element analysis and computation are significant. Hence, we want to study the convergence 
and complexity of AFEMs for multiple eigenvalue problems and focus on the following elliptic 
eigenvalue problems: find A G K and u £ H^{n) such that 



where A,c are coefficients stated precisely in Section [51 

We see that the analysis technologies for the convergence rate and complicity of AFEM in lit- 
erature are only valid for simple eigenvalues and their corresponding eigenfunctions, it can not 
be applied directly to multiple eigenvalue cases. The difficulty lies in that in context of multi- 
ple eigenvalue cases, it is not practicable to figure out the discreted eigenfunctions obtained over 
different meshes so as to approximate the same exact eigenfunctions. As a result, the standard 
technology of measuring the error of every eigenfunction does not work well any more, which re- 
sults in the difficulty when analyzing the reduction for error of the approximate eigenfunction over 
two consecutive meshes. Instead, we employ some measure between the eigenspace, the space of 
eigenfunctions, and its approximation. By using the similar perturbation argument in [TH [20] to- 
gether with eigenspace approximation technology and its crucial property that is also shown in this 
paper (see Lemma l3.ip . we obtain the convergence rate and quasi-optimal complexity of AFEM for 
multiple eigenvalue problems. 

The paper is organized as follows. In the next section, we shall describe some basic notation 
and review the existing results of finite element approximations to a class of linear second order 
elliptic boundary value and eigenvalue problems, which will be used in our analysis. In Section 
[31 we construct the a posteriori error estimators for finite element eigenvalue problems from the 
relationship between the elliptic eigenvalue approximation and the associated boundary value ap- 
proximation and then design adaptive finite element algorithm for the elliptic eigenvalue problems. 
We analyze the convergence and quasi-optimal complexity of the adaptive finite element eigenvalue 
computations in Sections [H and [H respectively. We present several numerical examples in Section 
[HI to illustrate our theory. Finally, we remark how our main results can be expected for the Steklov 
eigenvalue problems and inexact numerical solutions. 

2 Preliminaries 

Let C ]R'^((i > 1) be a polytopic bounded domain. We shall use the standard notation for Sobolev 
spaces T/F'''P(ri) and their associated norms and seminorms, see, e.g., [Din!- For p = 2, we denote 
H^ifl) = W^'^in) and H^ifl) ^ {v e H^{n) : v \dn= 0}, where v \oq= is understood in the 
sense of trace, || • \\s.n = || ■ ||s,2.si, and (•,•) is the standard inner product. Throughout this 
paper, we shall use C to denote a generic positive constant which may stand for different values at 
its different occurrences. For convenience, the symbol < will be used in this paper. The notation 
that A < B means that A < CB for some constant C that is independent of mesh parameters. All 
the constants involved are independent of mesh sizes. 

Let {Th} be a shape regular family of nested conforming meshes over fi: there exists a constant 
7* such that 



V • (AVu) +CU = \u in fi, 

||"||o,f2 = 1 on 
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where, for each T G Th, hx is the diameter of T, and pr is the diameter of the biggest bah contained 
in T,h = max{/iT : T G Th}- Let Eh denote the set of interior faces (edges or sides) of Th. 

Let S^'''^{fl) be a space of continuous functions on fl such that for v G S'^'''{^1), v restricted to 
each T is a polynomial of degree not greater than k, namely, 

S'^^'^ifl) ^{ve Cin) : v\t ePr VT G %}, 

where is the space of polynomials of degree not greater than a positive integer k. Set SQ'''{fl) — 
S'''''{n) n H^i^). We shaU denote S^^'^in) by S^iQ) for simplification of notation afterwards. 

2.1 A linear elliptic boundary value problem 

In this subsection, we shall present some basic properties of a second order elliptic boundary value 
problem and its finite element approximations. 

Consider the homogeneous boundary value problem: 

Lui = ft in Q, i ^ 1, - ■ ■ , N, 

(2.1) 

Ui = on d^l, 

where iV is a positive integer, L is a linear second order elliptic operator: 

Lu = - V • (AVu) + cu 

with A : ^ M.'^^'^ being piecewise Lipschitz over initial triangulation and symmetric positive 
definite with smallest eigenvalue uniformly bounded away from 0, and < c e L°°{n). 
The weak form of ((27T|) reads: find U = {ui, - ■■ , un) E (i?(J(ri))^ such that 

a{u,,v,) = b{f,,v,) yv,eH^in),i^l,--- ,N, (2.2) 

where 

a{u,v) ~ (AVujVv) + {cu,v) and b{u,v) = {u,v). 
We observe that a(-, •) is a bounded bilinear form over H^{U): 

\a{w,v)\ < Ca\\w\\i^n\\v\\i,n Vu;,u G i?o(ri) 
and for energy norm || • ||a,f2, which is defined by ||w||a,si = \^ a{w, w), there hold 

Ca||w||i,o < \\w\\a.n < Ca\\w\\i^n Vw e H^in), 
where Ca and Ca are positive constants. We understand that (|2.2p is uniquely solvable for any 

f,€H-^n){i = i,--- ,N). 

For i^(ri) with || • Ht fj ~ y/b{-, •), we see that there is a unique compact operator K : L'^{VL) ^■ 
7Jq(57) satisfying 

a{Kw,v) = b{w,v) yw e L^{n),v E H^i^). 
Define the Galerkin-projection Rh : Hq{Q) Vh = Sq{Q) by 

a{u - RhU, v) = Vu e H^{n) Vv G V^, (2.3) 

then apparently 

\\Rhu\\i,n<\\u\\i,n yueH^in). 
The following results can be found in [2j [34] . 
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Proposition 2.1. Let 



p,,{h) ^ sup mi \\Kf~v\\i^n. 

/eL2(0),||/||i,.n=l"e^'' 



Then (h) — > as /i — > and 

\\u - Rhu\\b,n < Pii{h)\\u - Vu G i?o(f^)- 

A standard finite element selieme for (12.21) reads: find Uh = (uui, • • • ,UN,h) S Vj^ satisfying 

a{u,^h,v^) =b{f,,v^) \fv,£Vh,i^l,---,N. (2.4) 

Let T denote tlie class of all conforming refinements by bisection of Thq ■ For 7^ G T, define the 
element residual TZxiui^h) and the jump residual Jsiui^h) by 

'R-riuiji) = fi- Lui^h = .fi + ^ ■ [AVui^h] ~ cui^h inTeTh, 
Jsiuiji) = -AVu'Ij^ ■ - AVtiJ,^ ■ = [[AVuiji]]£; • VE on E G Eh^ 

where E is the common side of elements and with unit outward normals and , 
respectively, and ue = i'~ ■ Let ujt be the union of elements sharing a side with T and be the 
union of elements which shares the side E^ that is, w^; = r+ H T~ . 
For T eTh,^e define the local error indicator f}h{ui^h^T) by 

fll{uiM,T) ^ h^\nT{uij,)\\lr^ + ^ ^Bps(ui./i)llo,£; 
and the oscillation osCh{ui,h, T) by 



OScl{Ui,h,T) = /l|||7^T(ui,/i) - iZT{uiM)\\l,T + X! ^sll^^^lwi,/!) - JE{uiM)\\l,E^ 

Ee£h,EcdT 

where w is the L^-projection of w G L'^(n) to polynomials of some degree on T or E. 
We define the error estimator fih{ui^h,^) and the oscillation osch{ui^h,^) by 

?7^(iiiji,ri) = ^ fil{ui^h,T) and osc^(wi,ft, r2) = ^ osc\{u^^h,T). 
TeTh,Tcn TeTh,Tcn 

For any [/^ = • • • , Miv,/i) e V^^ we set 

Af N 

VhiUh,^) ^^vKu^^h,^) and osc^([//i, rj) = ^ (^^^(m,,;,, fi). 

i=l 1=1 

In our analysis we need the following lemma [7]. 

Lemma 2.1. There exists a constant C* which depends on A, regularity constant j* , and coefficient 
c, such that 

osch{v,T) < osch{w,T) + C^,\\v - w\\i^^^, VwjWGF/i. (2.5) 
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We have the standard a posteriori error estimates for the finite element approximation of bound- 
ary value problems (PTTj) as follows (c.f., e.g., [HI UHl 133] ) 

\\ui - Ui^hWa^n < Cif]h{ui^h, (2.6) 

Clf|l{u^^h,^) - CloScliUi^h, f^) < \\Ut - M»,/i||a,n: (2-7) 

where Ci , C2 and C3 are positive constants depending on the shape regularity of the mesh 7^ . 
The adaptive algorithm for solving (12.41) without oscillation marking can be stated as follows. 





Algorithm Co 


Choose a parameter < 6 < 1. 


1. 


Pick any initial mesh T^q . 


2. 


Solve the system on Tha for the discrete solution Uhg = {ui,ho^ ' ' ' I'^N.ho)- 


3. 


Let fc = 0. 


4. 


Compute the local indicators fjh^. ■ 


5. 


Construct Aih^ C Th^, by Marking Strategy Eq and parameter 9. 


6. 


Refine 7/ij. to get a new conforming mesh Th^+i by Procedure REFINE. 


7. 


Solve the system on Th^+i to get the discrete solution Uh^+i = 




+ • ■ ■ ,unm + i)- 


8. 


Let k = k + 1 and go to Step 4. 



Marking Strategy Eq in Algorithm Co was introduced [131 [55] when = 1. It is used to 
enforce error reduction and can be defined as follows. 



Marking Strategy Eq 

Given a parameter < 6 < I. 

1. Construct a minimal subset of Th,. by selecting some elements in 
Tiif. such that 

E vUu>,,,T)>er)l{u>,,,n). 

2. Mark all the elements in Aih^- 



5 



The procedure REFINE here is not required to satisfy the Interior Node Property of [ITJ [25] • 
Given a fixed number b > 1, for any T/j.^ G T and a subset Mh^ C Th^. of marked elements, 

=REFINE(r,,,X„ J 

outputs a conforming triangulation Th^+i G T, where at least all elements of Mh^ are bisected 
6 times. We define RVh^^Th^^-^ = Th^XiThk ^ Thk+i) ^ the set of refined elements, thus Mhk C 

Lemma 2.2. Assume that Tho verifies condition (b) of Section 4 in 132]. Let {Thk\k>a be any 
sequence of refinements of Tho where Thk+i is generated from Th^ by Th^+i = R-EFINE(7hi , -^/it) 
with a subset Aih^ C Th^- Then 

k-l 

3=0 

is valid, where the hidden constant depends on Tho ^'^'^ ^■ 

By some primary operation, we can easily extend the results about the case of = 1 in ffj to 
vector problems, and obtain the following results. 

Theorem 2.1. Let {Uhk}keNo &6 sequence of finite element solutions of boundary problems pro- 
duced by Algorithm Cq. Then there exist constants 7 > and ^ G (0, 1), depending only on the 
shape regularity of meshes, the data, and the parameters used by Algorithm Co, such that for any 
two consecutive iterates k and k + 1 we have 

\\U- Uh,^, Win + iUh,^^ ' < {\\U - Uh, win + IVL iUh, , ^)) ■ 

Indeed, the constant 7 has the following form 

with some constant 5 G (0, 1). 

Lemma 2.3. Let Uh^^i G Vh^ and Uh^^^^^i G Vh^^-^ {I = 1, - ■ ■ , N) be discrete solutions of {2.4-^ over 
a conforming mesh Th^ and its refinement Th^^i with marked element Mh^- ^ = ^^''k^'^'^k+i 
be the set of refined elements, then the following localized upper bound is valid 

WUhk - Uhk+A\ln < Cl J2 vLiUhk,T), 
Ten 

where Uh^ = ••• ,Uh,^q) and Uhk+^ = ••• 

Proposition 2.2. Let Uh^.i G Vh^ and Uhk+i,i G Vh^+iil = 1, • ■ ■ ,N) be discrete solutions of {2.4^ 
over a conforming mesh Th^ its refinement Th^^i with marked element Mh^ ■ Suppose that they 
satisfy the energy decrease property 

WU- Uhk+, win + 7ooTc^h.+, {Uhk^^ , ^) < ^oiWU - Uh^ Win + 7oo^4 i^h^ , ^)) 

with 70 > being a constant and G (0, ^). Then the set TZ = Rj-f^^^j-f^^ ^ satisfies the Dorfler 
property 

Ten TeTh^ 
^^th = p^^^ff,';,'!;..)-^) . where Co = max(l, f ). 
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2.2 A linear eigenvalue problem 

A number A is called an eigenvalue of the form a(-, •) relative to the form b{-, •) if there is a nonzero 
function 7^ u G HQ{n), called an associated eigenfunction, satisfying 

a{u,v) ^ Xb{u,v) yveH^in). (2.9) 

We see that (|2.9p has a countable sequence of real eigenvalues 

< Ai < A2 < A3 < • • • 

and corresponding eigenfunctions 

Ul,U2,U3, ■■ ■ , 

which can be assumed to satisfy 

b{u,,Uj) = S^J, i,j = 1,2, • • • 

In the sequence {A^}, the A^ 's are repeated according to geometric multiplicity. 

The following property of eigenvalue and eigenfunction approximation is useful (see [51 135] ). 

Proposition 2.3. Let {X,u) be an eigenpair of S2.9\) . For any w € -ffo(i^) \ {0}, 

a{w,w) ^ a{w — u,w ~ u) ^b{w~u,w — u) 
b{w,w) b{w,w) b{w,w) 

A standard finite element scheme for (|2.9p is: find a pair of {Xh,Uh), where Xh is a number and 
^ Uh e Vh, satisfying 

a{uh, v) = Xhb{uh, v) yv e Vh. (2.10) 
Let us order the eigenvalues of (|2.10l) as follows 

< Ai,ft < >^2M < ■ ■ ■ < A„^,/i, Uh = dim Vh, 
and assume the corresponding eigenfunctions 

Ul^h,U2,h, ■ ■ ■ , Un,^,h 

satisfy 

b{u^,h,Ujji) ^S^j, i,j ^1,2,- ■■ ,nh. 

As a consequence of the minimum-maximum principle (see [5] or [S]) and Proposition 12. 3[ we 
have 

Xi < X^^h < K + C,\\u, - u^^hWl, i ^ 1,2,- ■■ ,nh. (2.11) 

Let A be any eigenvalue of (|2.9p with multiplicity q and M (A) denote the space of eigenfunctions 
corresponding to A, that is 

M(A) — {w £ ^o(^) : w is an eigenvector of (|2.9|) corresponding to A}. 

We see that A will be approximated from above by the Galerkin approximate eigenvalues: 

A < A/1,1 < • • • < Xh,q- 

Set 

(5/1 (A) = sup inf \\w-v\\a.n- 

weM{\),\\w\\b_n = i 

The following results are standard and can be found in literature (see, e.g., [2l|3l|8] or [35]). 
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Proposition 2.4. Let A be any eigenvalue of i2. 9\) with multiplicity q and Uh,i,--- ,Uh.q with 
\\uh,i\\b,n = ^fl = ^, ' ' ' iQ) be the Galerkin eigenf unctions corresponding to A/i.i,--- ,Xh.q, respec- 
tively. There hold 

\\u~Ehu\\b.n<Pa{h)\\u~Ehu\\a,n Vu G M(A), (2.12) 
\\uh,i - Euh,i\\b,n < p,-,{h)\\uhj - Euh,i\\aM, (2-13) 

where E denotes the orthogonal projection of the energy space onto M{X), Eh the orthogonal pro- 
jection onto span{uh,i, • • • , Uh^q}- 

For convenience, here and hereafter, we let M/i(A) — span{M/i^i, • • • ,Uh^q}- 

3 Adaptive finite element method 

Here and hereafter we assume that A is some eigenvalue of (|2.9|) with multiplicity q and u e M(A), 
and {\h.i,Uh.i) e M X V/j is the associated finite element eigenpairs of (|2.10p with = 1{1 = 

1, • • • ,q) that satisfy (I^TH)) and ^J^. 

Let A'' = a{EhU, Ehu), we conclude from (|2.1ip and (|2.12p that 

\\u - Ehu\\b,n + |A - A'^l < C(p„ (h) + \\u- Ehu\\a,n)\\u - Ehu\\a,n- (3.1) 

Note that (|2.9p and (|2.10p can be rewritten as 

u = XKu, Uh = XhRhKuh, 

where Rh is the Galerkin projection defined by (|2.3p . 

We see from the definition of Eh that for any u e M (A), there exist some constants {a/i,/(w)}'^i 

satisfying J2l=i {othj{u)f = 1 such that Euu = X^Li o:h,i{u)uh,i- Define = J2i=i ah,i{u)Xh.iKuh,h 
and we obtain from a direct calculation that 

EhU^Rhw''. (3.2) 

Theorem 3.1. Given u € M{X), let r{h) = Pa{h) + \\u - Ehu\\a,n- Then 

\\u - Ehu\\a,n = \\w''- Rhw''\\a,n + 0[r{h))\\u - Ehu\\a,n- (3.3) 

Proof. We obtain from the definition of that 

u-w^ = XKu - ahd{u)Xh,iKuh^i 
1=1 

= XK(u - Ehu) + XK(^ ah,i{u)uh,i) - ^ ah,i{u)Xh,iKuhj 

1=1 1=1 

q 

= XK{u - Ehu) + ^ ah,i{u){X - Xh,i)Kuh.i. 
1=1 

Since 

9 

A-A'' = ^(aM(w))'(A-A,.,,), 
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we have 



W'^ah.Mi'^- >^h,i)Kuh^i\\aSi < iy^a/i,K")(^~ ^h,i)\ 
1=1 1=1 

which together with the fact \\K{u — Ehu)\\a,n < ||u — Ehu\\b,n and (|3.ip leads to 

\\u-w''\\a,n < Cr{h)\\u - Ehu\\a,n- (3.4) 

Note that ([X^ imphes 

u - EfiU = w^- RhW^ +u-w^. 
Hence we obtain (|3.3p from (|3.4p . This completes the proof. □ 

3.1 A posteriori error estimators 

Following the element residual TZriuij^) and the jump residual JE{ui.h) for (|2.4p . we now define an 
element residual TZT{Ehu) and a jump residual JsiEhu) for (|2.10p as follows: 

TZxiEhu) = X'^'Eh.u + S/ ■ {A\/Ehu) ~ cEhU in T eTh, 
JsiEhu) = -AW{Ehu)+ -1^+ -AW{Ehu)- -i^- 
= [[AWEhu]]E ■ i^E on Ee£h, 



where E , z/+ and v are defined as those of section [2Jl 

For T (^Th , '^e define the local error indicator rjh{EhU,T) by 

nl{EhU,T) = hl\\nT{Ehu)\\lT+ E yE{Ehu)\\l,E 

EeSh.EcdT 

and the oscillation osch{EhU^T) by 

osc^ T) = hlWUTiEhu) - 'RT{Ehu)\\l^T + hE\\JE{Ehu) ^1^{E^\\1^e- 

EeSh^EcdT 

We define the error estimator riii{EiiU,Q) and the oscillation osCh{EhU,^l) by 

r]l{EhU,n)= ^ T]l{EhU,T) and osc^(i;,iW, il) = ^ osc\{EhU,T). 
TeTh,Tcn TGTh,Tcn 

For any = • • ■ , € V^;f, we let 

VhiUh, ^) =Y f]l{uhj,n) and oscl{Uh, f^) = ^ os4(u;i,i, f^) 

and for any EhU = {E^ui, • • • , E]-,Uq) e V^', we set 

q 9 
7?^ ^)^Y. 'll(EhUi,n) and oscliEhU, n) oscjXEhUi, Q). 

1=1 1=1 
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Given ho G (0, 1), define 

f{ho) = sup r{h). 
he{a,ho) 

Theorem 3.2. There exist constants Ci,C2 and C3, whieh only depend on the shape regularity 
constant 7*, coercivity constant Ca and continuity constant C'a of the bilinear form, such that 

\\u- Ehu\\a,n<Cnjh{EhU,n) (3.5) 



CWh{EhU,n) - CloscliEhU,Q) < \\u-Ehu\\l,, (3.6) 
provided ho 1. Consequently, 

\x-\"\<vl{EHU,n) 



f]l{EhU,n)~oscl{EhU,n) < |A-A''|. 
Proof. Recall that Lw'^ — J2l=i '^h.i{u)Xh^iKuh i. We obtain from (|2.6p and (|2.7p that 

Ww'^ -Rhw''\\a,n<Cir]h{.EhU,^) (3.7) 

and 

ClijliEhU, n) - CIoscI{Ehu, n) <\\w''^ Rhw'^wl^- (3-8) 



Combining (|3.2p . (I3.3p . (I3.7P with p.8|) . we complete the proof. In particular, we may choose the 
constants Ci, C2 and C3 satisfying 

C^Ciil + CHho)), C2 = C2{l-Cr{ho)), Cs = C^il - CHho))- (3.9) 

□ 

3.2 Adaptive algorithm 

Recall that the adaptive procedure consists of loops of the form 



Solve Estimate Mark Refine 



We assume that the solutions of the finite dimensional problems can be solved to any accuracy 
e0iciently0 The a posteriori error estimators are an essential part of the Estimate step. In the 
following discussion, we use r]h(Uh, ^) defined above as the a posteriori error estimator. 



^ In fact, we have ignored two important practical issues: the inexact solution of the resulting algebraic system 
and the numerical integration. We remark the discussion about the inexact solution in Section [71 
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Algorithm C 




Choose a parameter < ^? < 1. 




1. 


Pick an initial mesh Tho with mesh size /iq- 




2. 


Solve the system on Tho for the discrete solutions {Xho,U''J'ho, 
1, •••,?)• 


){l = 


3. 


Let fc = 0. 




4. 


Compute the local indicators r}h^ {uhk,i = 1, - ■ ■ ,q). 




5. 


Construct A^/i^ C Thk by Marking Strategy E and parameter 


9. 


6. 


Refine Thk ^ conforming mesh Th^+i by Procedure REFINE. 


7. 
8. 


Solve the system on Th^+i for the discrete solution {Xhk+i,i,Uhk+i 
Let k = k + 1 and go to Step 4. 


M = 



Our Algorithm C is nothing but Algorithm Co when Marking Strategy Eq is replaced by 
the following Marking Strategy E. 



Marking Strategy E 

Given a parameter < 6 <1. 

1. Construct a minimal subset Adh^ of Th^ by selecting some elements in 

such that 

E vlk{Uhk,T)>enliUhk,n). 

TeMuk 

2. Mark all the elements in Adh^- 



We shall now present the following property of eigenspace approximation that will play a crucial 
role in our analysis. 

Lemma 3.1. Given constant ^ e (0, 1). // 

J2 VH{UH,T)>eri%iUH,n), (3.10) 

TeMh 
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then for any orthonormal basis {ui}f^-^ of M{X), there holds 

^ r]jj{EHU,T)>e'vjiiEHU,n), (3.11) 

TeMn 

where ErU = {Ehui, ■■ ■ , EHUq) and 0' = ^ G (0, 1). 

Proof. Wc denote Ehui as vh,i- Since A'/^f (A) = spa,n{uH,i, • • • , U'H,q}, we have that there exists q 
constants /ShjU = 1, • " " ) ?) satisfying X)j=i {I^hj)"^ — 1 such, that 

Q 

VH,i =^l3H,jUH,j, l = l,---,q. 

We may analyze as follows 

q q 

1=1 j=l 
1=1 j=i 

3 = 1 1=1 

which together with the fact that < 1 yields 

r^jj{EHU,n)<2qrjjj{UH,n). (3.12) 

Wc sec from the definition of Eh that orthonormal basis of Mh{X), namely, 

Mh{X) = spanjw^f 1, • • • ,UH.q} = spa.n{vH,i, ■ ■ ■ ,VH,q}- So tliere exist q constants PhjU = 
1, • • • ,q) satisfying J2'j=i {i^HjT = ^ ^^^^ 

9 

Uh,1 =^$H,jVH,j, l = l,---,q- 
3=1 

Similarly, from the following estimation 

TeMh 1=1 TeMh 3=1 

1=1 j=l TeMh 

= 2EEfe)' E r), 

7=1 1=1 TeMh 



we obtain 



^ r,l{UH,T)<2q ^ ^1{EhU,T). (3.13) 

TeMn TeMn 
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Combing ([XTU]) . ([51^ and ((XT^ . we arrive at 



TeMf 



T£Mh 

which is nothing but ((3?TTt with 6*' = □ 

4 Convergence 

Following Theorem 13.11 we can establish some relationship between the two level approximations, 
which will be used in our analysis for both convergence and quasi-optimal complexity. 

Lemma 4.1. Let h, H E (0, ho), be any orthonormal basis of M{X), A^'' = a{EHUi, Ehui), 

w"-^ ^Y.l=iaHAui)^H.iKuH,i andW" = [w"'^--- Then 

\\U - EnU\\a,n = \\W ^ RhW"\\asi + O{f{ho)) i\\U ~ EhUWa.n + \\U- ShC/|U.o) , (4.1) 
oschiE^^U, n) = ^ChiRhW", n) + 0(7=(/io)) i\\U - E,Ma,n + \\U- EnUWa^n) , (4.2) 



VhiEhU, Q) = iih{RhW, n) + O{r{ho)) {\\U - EhUWa.n + \\U- EHU\\a,n) ■ (4.3) 

Proof. It is sufhcient to prove that for any u £ M{X), — '^'^^iah,i{u)Xh,iKuh^i, and — 
S?=i OiH,i{u)\H.iKuH.i, the following equalities hold, 

\\u~Ehu\\a,n = \\w" -Rhw"\\aS-i + 0{f{ho)){\\u-Ehu\\aM + \\u-EHu\\a,n). (4.4) 

osch{EhU,n) = 5^Ch{Rhw" ,n)^0{f{hi:,)){\\u- Ehu\\a,^i + \\u- EHu\\a,^i) (4.5) 

and 

m.{Ehu, n) = fih{Rhw",n) + O{f{ho)) {\\u - Ehu\\a,n + \\u - EHu\\a,n) ■ (4.6) 
First, we prove (|4.4[) . We see that 

WRhiw"" ~w")+w" -u\\a,a < Ww'' ~w"\\a,n + \\u~w"\\a.n 

< \\u-w"\\a,n + \\u-w''\\a,n, 

which together with p.4p leads to 

\\Rh{w'' - u.^) + w" - u\\a,n < f{ho){\\u - EHu\\a,n + \\u - Ehu\\a,n)- (4.7) 
Note that ([3?2l) implies 

u - Ehu ^w" - Rhw" + Rh{w" ~w^) + u- w", (4.8) 
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we get dm from (|i?71) . 

Next, we prove (|4.5p . We obtain from Lemma [2. II that 

SS'chiRhiw'' -w"),n)<\\Rh{w" -w'')\\aM, (4.9) 
which together with p.4p and (|4.7p yields 

5Fc,,(i?,,(w^ - u.''), ft) < f{ho)i\\u - EHu\\a,n + \\u- Ehu\\a,n)- (4.10) 
Due to EfiU — RhW^ + Rhiw^ — w^), we conckide from the definition of oscillation that 
^Ch{RhW^, n) = 5Ich{Rhw" + Rhiw'' - w"), n), 

then we arrive at (|4.5p . together with (I4.10p . 
We now prove (|4.6p . By (|2.7p . we have 

7y,.(i?,,(u;''-w;^),l]) < || (?«''- li;^) - i?h(u;'' ~ w;^)||a,o + 5Jc^(i?/,(u.'' - ^i) 
< h - w''\\a,n + \\u- w"\\a^n + \\Rhiw'' ~ w")\\a,n, 

where (|4.9p is used in the last inequality. Using p.4p and (j4.7p . we obtain 

77,,(i?,,(w''~w^),r!) < f(/io)(||u-^Hw||a,J2 + ||w--B,,u|U^n). (4.11) 

From (|4.1ip and the fact that 

jlh{Rhw'\ n) = 77/, (i?,,?/;^ + Rhiw'' ~w"),n), 

we get 

flh{RhW^\ 9) = fjhiRhw", n) + O{?{ho)) {\\u - Ehu\\a,n + \\u - EHu\\a,n) , 

which is nothing but (|4.6p since r]fi{RhW^, ft) — rih{EhU, fl). □ 

We are now in the position to present and analyze the error reduction result. 

Theorem 4.1. Let A e M 6e some eigenpair of i2.9\) with multiplicity q, {ui}^^-^ he an orthonormal 
basis of M{X), and {{Xh^j, Uh^^i), Z = 1, • • • , gjfcgNo « sequence of finite element solutions produced 
by Algorithm C of Section\^ Then there exist constants 7 > and a € (0, 1), depending only on 
the shape regularity of meshes, Ca and Ca, the parameter used by Algorithm C, such that for any 
two consecutive iterates k and k + \, we have 

\\U- En.^Mln + 7<+, {Eh.^.U, n) < a' {\\U - Ef^Mln + IvliE^^^U, n)) (4.12) 

provided ho <ti 1. Therefore, Algorithm C converges with a linear rate a, namely, the n-th iterate 
solution {X^'^'\Eh^ui){l — 1, - ■ ■ ,q) of Algorithm C satisfies 

\\U~E^,Mln + lvlAE^r.U,n) <Coa2" (4.13) 

and 

A'*"''-A<a2", (4.14) 
where Cq = ||C/ - E^Mln + IvliE^oU, n). 
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Proof. For convenience, we use {Xh,i,u,ij), {Xh,i,uh,i) to denote (Xh^^i^hUh^^-^j) and {Xhk,iiUh^,i), 
respectively. We see that it is sufficient to prove 

\\U - EhUWI^ + iigiEhU, n) < a'{\\U - EnUWln + Wh{EhU, n)). 

We derive from Lemma 13.11 that Marking Strategy E impfies that there exists a constant 
61' = ^ g (0, 1) such that 

J2 vhiEHU,T)>9'r^l{EHU,n). (4.15) 

TgMh 

Recall that w'*'' = X^Li C(hAui)^h.iKuhA, w"-^ = aH,t{'^i)XH,iKuH,i, we get from (|4.15p that 

for = {w^'^, • • • , w^''^), Marking Strategy Eq is satisfied with 9 = 6' . So, we conclude from 
Theorem 12. II that there exist constants 7 > and ^ e (0, 1) satisfying 

\\W" - RhW"\\l,, + jvliRhW", n) < e\\W" - EhUWI^ + ^HiEnU, n), (4.16) 

where 113. 2p is used when h is replaced by H. 

By Lemma [4. II and p.4p . we have that there exist constants Ci > 0, 6*2 > such that 

\\U-EHU\\ln+j7jjj{EHU,n) = \\W" - EHU\\l,, + ^7fH{EHU,n) 

+C,r^{ho)\\U -EhUWI^ (4.17) 

and 

\\U~EnU\\l,,+jriUEhU,n) = WW"" ~RhW''\\l^+jvliRHW",n) 

+C2f'{ho){\\U - E^UWl,, + EhUWI,,). (4.18) 
Combing (|4TT6ll . (|4T7l) with (|4?T8l) . we get when ft-o < 1 that 

\\u- EhUWI,, + / . vliEkU, n) 
, f(i -c.f (>.„)) + c.f(/.o) , _f ^,|,(E„t/,n). 

Since /iq <C 1 implies f(/io) ^ l, there exists 5i £ (0, 1) such that |(7i7^(/io)| < if ^0 ^ 1- Hence, 

\\u~EhU\\1,, + I . 4iEhU,n) 

1 - C27^(fto) 

when /iQ <C 1. Note that the constant a defined by 

\ l-C2i^{ho) ) 
satisfies a E (0, 1) if ho <C 1. We finally arrive at (|4.12p by using the fact that 



^2(l+5i)+C272(/i0) 

where 



< 7, 



7= ^- . (4.19) 

l-C2^Hho) 

This completes the proof. □ 
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5 Complexity 



We introduce a function approximation class as follows 

At,^{v e H : \v\s,^ < oo}, 

where 7 > is some constant, 

kk7 = supe inf {=f^T - ^Tho)' 

6>Q {rcno:inl{\\v-^VT\\la + h+i)osc?^ivT,T}V/^<e} 

and T C Thg means T is a refinement of Th^ . It is seen from the definition that, for all 7 > 0, 
— Al- For simplicity, here and hereafter, we use A^ to stand for Af, and use \v\s to denote 
So A'' is the class of functions that can be approximated within a given tolerance e by 
continuous piecewise polynomial functions over a partition T with number of degrees of freedom 

#r-#r^o <e-l/1^;|y^ 

To give the proof of optimal complexity of Algorithm C for the eigenvalue problem (|2.9p . 
we need some preparations. We know that in each mesh Th^, w'*'-'^' = 'Yl'i=i^hk.i{ui)^hi,.iKuh^,^i 
satisfies 

a{w'^'''\v) ^ (^ah^,i{ui)Xh^,iUh^,i,v) VweVh^, / = (5.1) 
i=i 

Thanks to Theorem l3.1l and Lemma [4.1l and their proofs, we are able to analyze the complexity of 
adaptive finite element method for multiple eigenvalue problems by using the complexity result for 
boundary value problems, which is similar to what was demonstrated in the convergence analysis. 

Using the similar procedure as in the proof of Theorem 14.11 we have 

Lemma 5.1. Let A G M &e some eigenvalue of i2. 9\} with multiplicity q, and {ui}J^^ be an orthonor- 
mal basis of M{X). Let {Xhk,iiUh^,i) G R x V/i^ and (Xh^+i.tiUh^+i.t) e M x Vh^^^ {I = 1, - ■ ■ ,q) 
be discrete solutions of \2.1U\) over a conforming mesh Th^ its refinement Th^+i with marked 
element A4hk- Suppose they satisfy the following property 

l2 I „„„2 

k+i 



\\U - Sh.+i c/lllo + l*osci^^, {Eh,^, u, n) 

where 7* > 0,/3* > are some constants. Then for the associated boundary value problem i5.1\) . 
we have 

WW''' - fi/..+iW^'1la,n + 7*oT4+,(^/..+,v^^'^^^) 



with 



>*'(l + <5i) + C'37'(/io)V^' 



\-C^^^(ha) / ^ l-C372(/^o) 



(5.2) 



where C5 is some constant depending on A , C and C*, i5i G (0, 1) is some constant as shown in 
the proof Theorem\^.l\ 
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Proof. By Lemma Wl] and (|3.4p . we have that there exists constant 6*3 > such that 

\\W" - RHW"\\l,,+^,5¥cl{RHW",n) = \\U-EHU\\l,,+j,5S'cl{EHU,n) 

+C^r\ho)\\W" -RhW" Wla, 

and 

\\W" -RhW"\\l^ + j,^cl{RhW",n) = \\U-EhU\\l^ + j,oscl{EhU,n) 

+C3P'{ho){\\W" - RhWWln + \\W"- RhW^WI^). 

Proceed the same procedure as in the proof of Theorem l4.11 we get 

\\W" - RhW"\\l,, + ^SS-cliRhW", n) < ~ RhW"\\1^ + %5S-cI{RhW", n)) 

with 



where C3 is some positive constant and i5i G (0, 1) is some constant as shown in the proof Theorem 
14.11 This completes the proof. □ 



The following statement is a direct consequence of E^^ui = Rh^w^'"^ , Proposition 12.21 and 
Lemma 15.11 

Corollary 5.1. Let A G M &e some eigenpair of V2. 9\) with multiplicity q, and {ui\f^^ he an 
orthonormal basis of M(X). Suppose that they satisfy the decrease property 

\\U - Eu,^Mln + l*oscl^^^{En,^,U,^) 

< l3l{\\U - EhMln+l*oscl{Eh,U,a} 



with constants 7* > and <E (0, y^). Then the set TZ = Rjj^^^j-f^^^^ satisfies the following 
inequality 

Ten "^sTfc^ 

with 9 ~ A ^^2M\- '^'^'^ ^0 = max(l, where /3* and 7* are defined in i5.2\) with 62 

being chosen such that (5^ G (0, i). 

Lemma 5.2. Let A G M &e some eigenvalue of 112. 9\) with multiplicity q, {ui}^^^ he an orthonormal 
basis of M{\), ui G A^{1 — I,-- - ,q) and Th^ he a conforming partition obtained from Tha- Let 
Thkj^i be a mesh created from Th^. upon making the set A4hk according with Marking Strategy E 

with 9 G (0, c'\C^+{WCic'^)i) ^- {{^hkd,Uh^,i)-,l = 1, • • • ,<7}fcGNo be discrete solutions of \2.10\l 
over a conforming mesh Th^ o,nd Mhf.{X) — span{uh^,,i, • • ■ ,Uht,q}. Then 



#M^,,<c(\\U^E,,Mln + loscliEn.U,n)) (9*-^ ^ l^i/l^^) , (5.3) 

1=1 

where the constant C depends on the discrepancy between 9 and ^■A^ij'i^^-^^2c'^c'^)-y) ■ 
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Proof. Let ^,^1 e (0, 1) satisfy /3i S (0,/3) and 

''^ C|(C2 + (1 + 2C2C2)7)^' ^ 

Choose 

e = (11 [/ ~ S,, (7|12 + 7os4^ (7, 17)) 

and let T/i^ be a refinement of 7ho with minimal degrees of freedom satisfying 

- £;,j^M,||^f2 + (7 + l)os4j£;,i^U;,17) < — , / = I,-- - 

which means 

r - E,,Mls-i + (7 + l)osc2^(ii;„,C/, 1]) < e\ (5.4) 
We get from the definition of that 

which implies 

#7/1^ - #77io 



(=1 



Let Thk + be the smallest common refinement of Tk^ ^-nd Th^ ■ Note that both Th^ and 7^^ are 
refinements of T^^, we have that the number of elements in T^t + that are not in Th^ is less than 
the number of elements that must be added to go from 7^^ to Th^ , namely, 

Let w^^''- = J2i=i oth,,i{ui)\h,,iKuh,,i = K[Yh=i aK,i{'^i)>'K,iUh,,i) = 1, • • -g), namely 

9 

We obtain from Lemma EH and the Young inequality that 

< 25Icl^{Rh,W''%n)^2Cl\\Rh,,^W''^ - RKW^'Wln- 

Note that 7^^ ^ is a refinement of 7^^ , L^-projection error are monotone and the following 
orthogonality 

is valid, we arrive at 

Ww'^^ ^ R^^^w''^\\l^ + ^5S-cl^ jR^^^w''^ ,n) 

< WW''^ - R,,^W''^\\1,, + ±oscl{Rh^W''%n). 
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Since (|2.8p implies 7 < we obtain that a = ~ j £ {0,1) and 

< WW'^^ -R^^w'^^\\l,, + ±o.scliRH,W'^^,n) 

< WW'^^ -RH^W'^^\\l,, + i^ + a)oscl{R^^W'^^,n). 
Applying the similar argument in the proof of Theorem 14.11 we may conclude that 

< f3'o{\\U- E,,Mln + (7 + o)oscl{En^U, 17)) 

< Pl{\\U-EhMln + {l + ^)oscl^{EKU,^)), (5.5) 

where 

(\ 1/2 
l + 5^ + C2i\ho) \ 

and Si is the constant appearing in the proof of Theorem 14. II Combining (|5.4p and (|5.5p . we then 
arrive at 

\\U- + loscl^ jE^,,^U,n) <P'{\\U~ Eh,U\\l+loscl{E^,U,n)) 

with $ = ^/3o/3i- 

Let (5i G (0, 1) be some constant satisfying 



(1 + <^i)'/3? < /3^ (5.6) 

which implies 

(l + 5i)/??<l. (5.7) 
We see from ho 1 and (15. 7p that e (0, 5). Thus we get from Corollary 15.11 that Th^ + satisfies 

^6^'^..+ Ten, 

^ ^ / /3^(1 + ^i) + C'37^(M \ 

1^ l-C372(/jo) j 

From the definition of 7 (see (14.191) ) and 7 (see (|2.8p ). we obtain that 7 < 1. On the other hand, 

- (=,2 _ _ 

we have C3 > 1 and hence Cq ~ Consequently, we can write 6 as 9 — 



Cj(l-2/3") 



Ci(^ + (1+2CJC?)) 

Since /iq ^ 1, we obtain that 7 > 7 and /3 G (0, from (|5.6II . Using p.9p . we get that 



> —(1-/3^) 



C|(^ + (1 + 2C^,Cl)) Cii^ + (1 + 2C2C?)) 



c 



(l-CT(.o))^ .(1-/32), 



(1-Cf (fto))^ Hl+Cf (/io))27 ^ * {l + Cf{ho)) 
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which together with the fact that ho <^ 1 and 7 > 7 yields 

- c|(^ + (1 + 2C^Cl)) ^1(^1 + (1 + 2C2C2)7) 

Therefore, we have 

#M,,, < #7^ < #r^,,^ - < - #%o 



This is the desired estiraate (j5.3p with an exphcit dependence on the discrepancy between 9 and 
c'^(C'^+{i+2C^c'^)i) • This completes the proof. □ 

We are now ready to show that Algorithm C possesses quasi-optimal complexity. 

Theorem 5.1. Let A G M 6e some eigenvalue of i2. 9\) with multiplicity q, {ui}1^^ be an orthonormal 
basis of M{X), and ui G A'^(l = 1, • • • , g). Let {{Xh^i,Uh^..i), I = 1, ■ ■ ■ , gjfcgNo ^6 ^ sequence of finite 
element solutions produced by Algorithm C of Section\3l and Mii,{X) = span{ufi,^i, ■ ■ ■ ,Ufi^ q}. 
Then the n-th iterate solution space Mh^{X) 0/ Algorithm C satisfies the quasi-optimal bound 

\\U- E,,^U\\l^ + joscljE,,^U, n) < (#77.„ - Wh„r'% 
where the hidden constant depends on the exact solution (A, M(A)) and the discrepancy between 9 



and 



Proof. We see from p. 61) that 

\\U- E,,Mln + IvliEh.U, n)<C{\\U- E,Mln + loscl^En.U, f!)), 
where C — max(l + ^). Hence, we get from ()5.3p that 

Note that Theorem 14.11 implies 

\\U ~ E^,^Mln+lvL^AEk,+.U,n) < a'{\\U-E^Mln + Jvl{Eh,UM))- 
Thus for < fc < n, we arrive at 

-l/2s 

-l/2s 



{\\U - EnMln + lvU^h^U^n) 
< a("-^)/^ (lie/ - E^^UWl^ + 7< {Eh^u, Q)) 
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We next employ Lemma [2.21 to deduce that 

-l/2s 



ri-1 



1=1 k=0 



— 1/2 " 



1=1 k=l 



which together with the fact a < 1 leads to 

-l/2s - « 



1=1 

Since osc{Eh^U,n) < rjh^{Eh^U,ri), we conclude 



This completes the proof. □ 

6 Numerical examples 

In this section, we show some numerical experiments for both linear finite elements and quadratic 
finite elements in three dimensions to illustrate the theoretical results obtained in this paper. 

Our numerical experiments were carried out on LSSC-III in the State Key Laboratory of Sci- 
entific and Engineering Computing, Chinese Academy of Sciences, and our codes were based on 
the toolbox PHG of the State Key Laboratory of Scientific and Engineering Computing, Chinese 
Academy of Sciences. 

For the convenience of analysis for our numerical results below, we define the following notations. 
Let A € R be some exact eigenvalue with multiplicity q, M{\) be the corresponding eigenfunction 
space as defined before, and {ui]1^^ be some orthonormal basis of M(A). Let {Xh,i,Uh,i)(l = 
1, • • • , g) be the finite element approximation, Mh{X) — span{uft,i, • ■ ■ , Uh,q}. For M{X) and Mh{X), 
we define the a posteriori error estimators rjh{Mh{X)^ fi) and the error D{M{X), Mh{X)) over mesh 
Th as follows: 

9 

7jhiMi,ix),n) = J27jhiuh,i,n), (6.1) 

1=1 

q 

D{M{X),Mh{X)) = Y.W^'i - ^hUiWa.n. (6.2) 



We see from the proof of Lemma 13.11 that 

9 

ilh{Mh{X),n) ^Y.'1h{EhUun). 



1=1 



Thus we get from Theorem 13.21 that 

i^h{Mh{X),n) ^ D{M{X),Mh{\)). 
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Example 1 Consider the following harmonic oscillator equation, which is a simple model in 
quantum mechanics |18j : 

-^Au+^\xfu = Xu mR^, (6.3) 



where |a;| = -^/jxip + |a;2p + la^sP- The eigenvalues of ()6.3p are A„ = n + ^ with multiplicity 
ra(n+i) _ 1,2,---), and its associated eigenfunction is m„ — 7e~l^l /^Hn{x) with any nonzero 
constant 7 and i?„(a;) = (-l)"e^^^(e~^^). 

Since the solution of (jG.Sp exponentially decays, we may solve it over some bounded domain fl. 
In the computation, we solve the following eigenvalue problem: find (A,m) e M x Hq{^1) such that 



\u 

n 



^dx ~ 1 and 



[ -iAu+i|a;pii ^ Xu in fi, 

I u = on dfl, 

where ft = (—5.5, 5.5)'^. We calculate the approximation of the first two smallest eigenvalues Ai and 
A2 with multiplicity 1 and 3, respectively, and their corresponding eigenfunction spaces M(Ai) and 
Af (A2) with dimension 1 and 3, respectively. 

Some cross-sections of the adaptively refined mesh constructed by the Marking Strategy E 
are displayed in Fig. 16.11 from which we observe that the mesh is denser in the center of the domain 
where the solutions oscillate quickly than in the domain far away from the center where the solution 
is smoother. This shows that our adaptively refined mesh can catch the oscillation of the solution 
efficiently and the a posteriori error estimators we designed are efficient. Our numerical results are 
presented in Fig. 16.21 and Fig. 16.31 Since the multiplicity of the first two smallest eigenvalues are 
1 and 3, respectively, for the discrete problem, we calculate the first 4 eigenpairs. We see from 
the left figure of Figure 16.21 that the convergence curves of error for all eigenvalues by using linear 
finite elements are parallel to the line with slope — |. Besides, we also observe that the convergence 
curves for the second, the third and the forth eigenvalues overlap together, this coincide with the 
fact that the multiplicity of the second eigenvalue is 3. Meanwhile, from the left figure of Fig. 
16.31 we see that by using linear finite elements, the convergence curves of the a posteriori error 
estimators for eigenfunction space rjh{Mh{Xi)) and rih{Mh{Xi)) are parallel to the line with slope 
— ^. Since 77?i(M/i(A), 51) « D{M{X)^ Mh{X)), we obtain that the convergence curves of error for 
the eigenfunction space D{M{Xi), Mh{Xi)) and -D(M(A2), Mh{X2)) are also parallel to the line with 
slope — ^ . This means that the approximation of eigenvalues as well as the eigenfunction space have 
reached the optimal convergence rate, which coincides with our theory in Section (5] We have the 
similar conclusion for the quadratic finite elements from the right figures of Fig. 16.21 and Fig. 

Example 2 Consider the Schrodinger equation for hydrogen atoms: 

1 ^ 1 

with / \u\'^dx — 1. The eigenvalues of (|6.4p are A„ = —^^{n = 1, 2, • • • ) and the multiplicity of 
Xn is v? (see, e.g., [I9]). 

Since the eigenfunctions of (16.41) decay exponentially, instead of (|6.4p . we may solve the following 

eigenvalue problem: find (A,m) e M x i/o(f^) such that / jupfia; = 1 and 

Jn 



- -5- ) -u = Am in (6.4) 



"^^"r)" ^ (6.5) 
u = on 9fi, 
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(a) Linear finite elements (b) Quadratic finite elements 

Figure 6.1: The cross-sections of an adaptive mesh of Example 1 
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(a) Linear finite elements (b) Quadratic finite elements 

Figure 6.2: The convergence curve of relative error for the eigenvalues of Example 1 



where J7 is some bounded domain in IR^. In our computation, we choose fl = (—20.0,20.0)^ 
and find the first 2 smallest eigenvalue approximations and their associated eigenfunction space 
approximations. Since the multiplicity of the n-th smallest eigenvalue is n^, for the discrete problem 
of (|6.5p . we calculate the first 5 smallest eigenvalues and their associated eigenfunctions. 

Fig. 16.41 is the cross-sections of the adaptively refined mesh constructed by Marking Strategy 
E. Similarly, we see that for both the linear finite elements and quadratic finite elements, the mesh 
is much denser in the center of the domain where the solution oscillates quickly than in the domain 
far away from the center where the solution is smooth. This means that the a posteriori error 
estimators we used are efficient. 

The numerical results are presented in Fig. 16.51 and Fig. 16.61 Similar to Example 1, Fig. 
16.51 shows that the convergence curve for all the eigenvalues obtained by linear finite elements and 
quadratic finite elements are parallel to the line with slope — | and — |, respectively, which means 
all the eigenvalue approximations reach the optimal convergence rate for both linear finite element 
and quadratic finite element. Meanwhile, we see from Fig. 16.61 that convergence curve for the a 
posteriori error estimators for eigenfunction space rih{Mh{Xi)) and rjh{Mh{X2)) obtained by linear 
finite element are parallel to the line with slope — \ , and those obtained by quadratic finite element 
are parallel to the line with slope — |. We observe that the approximation of eigenfunction space 
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(a) Linear finite elements (b) Quadratic finite elements 

Figure 6.3: The convergence curves of rjh{Mh{Xi), ^) and ■qh{Mh{\2),^) for Example 1 



(a) Linear finite elements (b) Quadratic finite elements 

Figure 6.4: The cross-section of an adaptive mesh of Example 2 



has also reached optimal convergent rate. These results validate our theoretical results. 

Example 3 Consider the following eigenvalue problem which is defined in a non-convex domain: 

find (A,m) e R X Hl{^) such that / \u\'^dx = 1 and 



-Au = 
2 



Au in n, 

on on, 



(6.6) 



where ft — (—5.0,5.0)^ \ (0,5.0)^, see the left figure of Fig. 16.71 below. We observe from the 
numerical calculation that Ai = 0.210651 with multiplicity 1 and A2 = 0.331779 with multiplicity 
2. 

The surface of the adaptively refined meshes constructed by Marking Strategy E is shown in 
Fig. 16.71 Besides, some cross-sections are displayed in Fig. 16.71 We see from these figures that for 
both linear finite elements and quadratic finite elements, the mesh is much denser along the lines 
where the solution is singular than in the domain far away from the singular lines. It indicates that 
our error estimator and marking strategy are efficient. 

Our numerical results are listed in Fig. 16.81 and Fig. 16.91 Similar to Example 1 and Example 
2, we can also see that the approximations of eigenvalue as well as eigenfunction have reached 
optimal convergence rate, which coincides with our theory in Section [51 
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10000 

■e of freedom 



(a) Linear finite elements 



10000 

e of freedom 



(b) Quadratic finite elements 



Figure 6.5: The convergence curves of relative error for eigenvalues of Example 2 




Figure 6.6: The convergence curves of 77/i(M^(Ai), fi) and rih{Mh{\2) , ^) for Example 2 



7 Concluding remarks 

We have studied the convergence and complexity for elliptic multiple eigenvalue problems. Now we 
turn to address how to apply the same arguments to the Steklov problem that consists in finding 
A e M and u ^ such that 

f -V • (AVu) +CU =0 in fl, 
1 (AVw) • it = Xuv on dfl, 

where it is the outward unit normal vector of fl on dfl. 

We set H ^ H\n), Vh ^ S''{n), W = L'^idfl), •) = (•, Oao and consider the non- 
homogeneous Neumann problem as a model problem as follows: 

Lui = in n, i = 1, • • • ,iV, 

(7.1) 

(AVu,) - it = on on. 
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(a) nonconvcx domain Q 



(b) Linear finite elements 



(c) Quadratic finite elements 



Figure 6.7: The surface of an adaptive mesh of Example 3 
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the first eigenvalue 
the second eigenvalue 
th^ third eigenvalue 
i^ith slope -4/3 




10000 

e of freedom 



10000 

e of freedom 



(a) Linear finite elements (b) Quadratic finite elements 

Figure 6.8: The convergence curve of relative error for eigenvalues of Example 3 



Define the element residual TZriuiji) and the jump residual JE{ui.h) for (|7.1I) as follows: 

JE{uiJi) 



j [[A\7u,^h]]E -^E onE e £h 
I ylVui^ft • it - /i on -Be 5/1, 



where Sh denote the set of boundary faces. For T e 7/i, we denote the local error indicator 
f}h{ui.h,T) by 



and the oscillation osch(ui,h,T) by 



X! hE\\JE{UtJi)\\l^E 



EeSh^Sh-EtZdT 



We see that Lemma [2.51 is also valid for (j7.2p . 



hE\\JE{uiji) - Je{uiJi)\\1^e- (7-2) 



Ee£hUSh,EcdT 
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(a) Linear finite elements (b) Quadratic finite elements 

Figure 6.9: The convergence curves of 77/i(Mft(Ai), fi) and ■qh{Mh{X2),^) for Example 3 



In context of Stekfov eigenvalue problems, we define 

nriEhu) = V ■{AVEhu)~cEhU in TeTh, 

j [[AVEhu]]E -VE on EeSh 

I AVEhU ■ it - X'^EhU onE e Sh, 



JE{Ehu) 



where n denotes the outward unit normal vector on E E Sh- For T £ Th, we define the local error 
indicator rih{EhU,T) by 



rii{EhU,T)^hl\\nT{Enu)\\lr- 
and the oscillation osch{EhU,T) by 



E 



^ \\JE{Ey,u)\\ 

%,E 



hE\\JEiEhu) - JE{Ehu)\\lE- 



oscl{EhU,T) = h^TlTiEhu) -TZTiEhu)\\lT + 

Ee£hUSh,EcdT 

We obtain by using the same argument that Theorem 14.11 and Theorem EH] are valid for the Steklov 
problem with multiple eigenvalues. 

In our numerical analysis above, for convenience, we assume that the algebraic eigenvalue prob- 
lem is exactly solved and the numerical integration is exact. Indeed, the same conclusion can be 
expected if all the numerical errors are taken into account, including both the error resulting from 
the inexact solving of the algebraic eigenvalue problem and the error coming from the inexact nu- 
merical integration. Suppose (A, u) is an eigenpair with the multiplicity of A being q, the exact 
solution on mesh Tk are {{Xh,i,Uh,i)}1^i, and the the solution considering the numerical error are 
{{Xh,i-,Uh,i)}1^i- If the numerical errors resulting from the solution of algebraic system and the 
numerical integration are small enough, say, satisfy 

q q 



i=l 



i=l 



with r(ho) <C 1 for ho ^ 1, then we have from the following triangle inequality 

\\Ui - Uh.i\\a < \\Ut - Uh.i\\a + \\uh,i - Mft.iHa, 

that our main results obtained in this paper hold true for inexact algebraic solution and inexact 
numerical integration, too. 
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